IL NUOVO CIMENTO 



Vol. ?, N. ? 



? 



An alternative derivation of the Lense-Thirring 
drag on the orbit of a test body 

Lorenzo Iorio( 1 ) 

( ) Dipartimento Interateneo di Fisica dell' Universitd di Bari, Via Amendola 173, 70126, Bari 



Summary. — In the weak field and slow motion approximation of general relativity 
a new approach in deriving the secular Lense-Thirring effect on the orbital elements 
of a test body in the external field of different central rotating sources exhibiting 
axial symmetry is presented. The approach adopted in the present work, in the 
case of a perfectly spherical source, leads to the well known Lense-Thirring formulas 
for all the Keplerian orbital elements of the freely falling particle. The corrections 
induced in the case of a central nonspherical, axisymmetric body are also worked 
out. 

PACS 95.30 - Relativity and gravitation. 

PACS 04.20 - Fundamental problems and general formalism. 



1. Introduction. 

Latest years have seen increasing efforts devoted to the measurement of the general 
relativistic Lense-Thirring effect [1, 2] in the weak gravitational held of the Earth by 
means of artificial satellites. At present, there are two main proposals which point to- 
wards the implementation of this goal: the Gravity Probe-B mission [3], and the approach 
proposed in [4] which consists in using the actually orbiting LAGEOS laser-ranged satel- 
lite and launching another satellite of LAGEOS type, the LARES, with the same orbital 
parameters of LAGEOS except for the inclination which should be supplementary with 
respect to it. At present, both these two satellites have not yet been launched: however, 
while the GP-B is scheduled to fly in 2002, to date the LARES mission has not been 
approved by any Spatial Agency. 

Recently Ciufolini [5] has put forward an alternative strategy based on the utilization 
of the already existing LAGEOS and LAGEOS II which could allow the detection of the 
Lense-Thirring drag at a precision level of 20 % [6, 7, 8]. While the GP-B mission is 
focused on the gravitomagnetic dragging of the spin of a freely falling body, in the LA- 
GEOS experiment it is the whole orbit of the satellite which is considered to undergo the 
secular Lense-Thirring precession. More exactly, among the Keplerian orbital elements 



© Societa Italiana di Fisica 



1 



2 



LORENZO IORIO 



the node and the perigee are affected by the gravitomagnetic force. For LAGEOS and 
LAGEOS II its effect consists in secular precessions amounting to: 



(1) ^LAGEOS „ 31 mas/y; 

(2) f^AGEOSII „ 31 5 mag/y; 

(3) ^ageosii „ _ 57 mag/y; 



where mas/y stands for milliarcseconds per year( 1 ) 

In this work it is presented an alternative strategy in order to derive the secular 
gravitomagnetic precessions which reveals itself useful especially in the prediction of the 
behavior of all the Keplerian orbital elements of the test body in the gravitational field 
of different kinds of central rotating sources. Indeed, the calculations involve not only 
the secular effects for a perfectly spherical source, which lead to eqs.(l)-(3), but also for 
a central body which is only axially symmetric. 

Until now the general relativistic motion in the field of a non spherical central body 
has been treated in [9, 10] for a static, nonrotating source and in [11, 12, 13] for a rotating 
source in the context of the GP-B mission. The influence of the non sphericity of the 
central source on the gravitomagnetic clock effect [14] is investigated in [15]. 

Our calculations are based on the use of the Lagrangian planetary equations [16] and a 
noncentral hamiltonian term whose existence, if from one hand can be rigorously deduced, 
from the other hand can be intuitively guessed by analogy from the corresponding term 
in electrodynamics for a charged particle acted upon by electric and magnetic fields. 



2. The gravitomagnetic potential 

In general, it can be proved that the general relativistic equations of motion of a test 
particle of mass m freely falling in a stationary gravitational field, in the weak field and 
slow motion approximation, are given by [15]: 



( J ; - m(E g + x B g ). 



In eq.(4) c is the speed of light in vacuo while E g and B g are the gravitoelectric and the 
gravitomagnetic fields, respectively. 

If a perfectly spherically symmetric rotating body is assumed as gravitational source, 
in eq.(4) E g = —GMr/r 3 is the Newtonian gravitational field of a spherical body, with 
M its mass and G the Newtonian gravitational constant, while B g is given by: 

^ , n G r J-3(J-r)r n 
(5) B 9 = VxA 9 .2-[ 

in which: 

G J x r 



(6) A 9 (r) ~ -2 



c r 6 



( x ) The perigee of the LAGEOS has not been considered since its rate is hard to measure [5]. 
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In eqs.(5)-(6) J is the angular momentum of the central body. The field A g , named 
gravitomagnetic potential, is due to the off-diagonal components of the spacetime metric 
tensor: 

(7) g^u = V + V" M> v = 0, 1, 2, 3, 

where r]^ is the Minkowski metric tensor and (A fl )fc = h k, k = 1,2,3. In obtaining 
eq.(6) a non rotating reference frame K{x, y, z} with the z axis directed along J and 
the {x, y} plane coinciding with the equatorial plane of the gravitational source has been 
used. The origin is located at the center of mass of the central body. 

The gravitomagnetic potential generates a non central gravitational contribute due 
uniquely to the angular momentum of the gravitational source that the Newtonian me- 
chanics does not predict, though the conditions of validity of eq.(4) are the same for 
which the latter holds as well. So it is possible to speak of mass-energy currents whose 
motion exerts a non central gravitational force on a test massive body analogous to the 
Lorentz force felt by a charged particle of charge q and mass m when it is moving in a 
electromagnetic field. Indeed, its equations of motions: 

(8) m^f ~ g(E + -v x B) 

are formally analogous to eq.(4). Eq.(8) can be derived by means of the Lagrangian: 

(9) C e . m . = \mv 2 -qV + ^-{wA), 

£ c 

where v is the velocity of the particle while V and A are the scalar and vector potential, 
respectively, of the electromagnetic field. 

Since one of the most promising way to detect the gravitomagnetic precession consists 
in employing artificial Earth satellites, it would be helpful to derive the rate equations 
for the change in the parameters that characterize their orbits. To this aim one could 
introduce "by hand" a perturbative term k (v • A fl ) in the gravitational Lagrangian of 
the particle and use it in some particular perturbative scheme; the constant k would be 
determined by means of dimensional considerations and taking in account that it should 
be built up of universal constants. In fact it is possible to show that a non central term 
analogous to | (v • A) can be rigorously deduced in the Lagrangian of a test body in the 
gravitational field of a spinning mass-energy distribution, and that it can be exploited in 
deriving straightforwardly the effect of the gravitomagnetic potential on the Keplerian 
orbital elements of the test body. 

3. The rate equations for the Keplerian orbital elements 

The relativistic Lagrangian for a free particle of mass m in a gravitational field can 
be cast into the form: 

(10) C = C (Q) +£ {1 K 

In eq.(10) the term D- 1 ' refers to the off-diagonal terms of the metric: 



(11) 



£(D = ™g ok x°xK 
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In this case, recalling that the slow motion approximation is used, eq.(ll) becomes [15]: 
(12) = C gm ~™ (A 9 -v). 

In this paper it is proposed to adopt C gm given by eq.(12), with A g given by eq.(6), in 
order to derive the Lense-Thirring effect on the orbital elements of a particle in the field 
of a rotating gravitational source. 

To this aim it must be assumed that under the gravitomagnetic force the departures of 
the test body' s trajectory from the unperturbed Keplerian ellipse are very small in time. 
This allows to introduce the concept of osculating ellipse. So the perturbed motion can 
be described in terms of unperturbed Keplerian elements varying in time. Consider the 
frame K{x, y, z} previously defined and a frame K {X, Y, Z} with the Z axis directed 
along the orbital angular momentum 1 of the test body, the plane {X, Y} coinciding 
with the orbital plane of the test particle and the X axis directed toward the pericenter. 
K{x, y, z) and K {X, Y, Z} have the same origin located in the center of mass of the 
central body. The Keplerian orbital elements are: 

• a, e 

The semimajor axis and the ellipticity which define the size and the shape of the orbit 
in its plane. 

• ft, i 

The longitude of the ascending node and the inclination which fix the orientation of the 
orbit in the space, i.e. of K {X, Y, Z} with respect to K{x, y, z}. The longitude of 
the ascending node Cl is the angle in the equatorial plane of K{x, y, z} between the x 
axis and the line of nodes in which the orbital plane intersects the equatorial plane. The 
inclination i is the angle between z and Z axis. 

• w, M 

The argument of pericenter and the mean anomaly. The argument of pericenter to is 
the angle in the orbital plane between the line of nodes and the X axis; it defines the 
orientation of the orbit in its plane. The mean anomaly M. specifies the motion of the 
test particle on its orbit. It is related to the mean motion n = 

(GM) 1 /2 a -3/2 thr0U gh 

M. = n(t — t p ) in which t p is the time of pericenter passage. 
It is customary to define also the longitude of pericenter 

• w = fi + UJ, 

the argument of latitude 

• u = w + f 

where / is the angle, called true anomaly, which in the orbital plane determines the 
position of the test particle with respect to the pericenter, and the mean longitude at 
the epoch to 

• e = ZD + n(to — t p ). If to = 0, it is customary to write e as Lq = vu — nt p . 

The matrix R x x for the change of coordinates from K {X, Y, Z} to K{x, y, z} is: 

(cos f2 cos to — sin f2 cos i sin u> — cos fi sin w — sin £1 cos i cos u 
sin fl cos uj + cos ft cos i sin uj — sin sin u> + cos Q cos i cos uj 
sinisinw sin i cos u) 



sin Cl sin i 
- cos Q sin i 
cosi 



Using eq. (13) andX = rcos/, Y = rsin/, Z = it is possible to express the geocentric 
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rectangular Cartesian coordinates of the orbiter in terms of its Keplerian elements: 

x = r(cos u cos f2 — sin u cos i sin fi) 

(14) ^ y = r(coswsin£l + sin u cos i cos fi) 

z = r sin u sin i. 

Redefining suitably the origin of the angle f2 so that cosfi = 1, sinSl = 0, the previous 
equations become: 

x = r cos u 

(15) ^ y = r sin u cos i 

z = r sin u sin i. 

Considering for the test particle the total mechanical energy with the sign reversed, 
according to [16], T = —Etot = —{T + U), where T and U are the kinetic and potential 
energies per unit mass, it is possible to work out the analytical expressions for the rate of 
changes of a, e, i, Q, u>, M. due to any non central gravitational contribution. To this 
aim it is useful isolating in hi the central part — C of the gravitational field from the terms 
—1Z which may cause the Keplerian orbital elements to change in time: U = —C — 1Z. In 
this way T becomes: 

-r- GM m ^ GM m 

16 T= +K-T=—- +K. 

r 2a 

Concerning the perturbative scheme to be employed, the well known Lagrange planetary 
equations are adopted. At first order, they are: 



(17) 
(18) 



- — — 

dt na dM' 

de _ 1 - e 2 8TZ (1 - e 2 ) 1 /2 QTl 



dt na 2 e 8A4 na 2 e du ' 

di _ 1 dll 1 dn 

( ' dt'^'na^l-e^sim IkJ ~ na 2 (l - e 2 ) 1 / 2 sini W 

dn i on 



(20) 



dt na 2 (l — e 2 ) 1 / 2 sini di 



dw _ _ . 1 dn (1-e 2 ) 1 / 2 dn 

dt na 2 (l — e 2 ) 1 / 2 sini di na 2 e de ' 

dM _ _ 1-e 2 dn 2 dn 
dt na 2 e de na da ' 

The idea of this work consists in using C gm to obtain a suitable non central term n gm to 
be employed in these equations. This can be done considering the Hamiltonian for the 
test particle: 

(23) H = p • v - C. 
Inserting eq.(10) in eq.(23) one has: 

(24) H = H {0) + H gm , 
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with H gm = — ^ (A g • v). So it can be posed: 

(25) ft ffm = _Z^i = I ( A 9 .v). 

Now it is useful to express eq.(25) in terms of the Keplerian elements. Referring to 
eq.(6), eq.(15), recalling that in the frame K{x, y, z} J = (0, 0, J) and that for an 
unperturbed Keplerian motion: 

r a(l - e z ) 

it is possible to obtain, for a perfectly spherical central body: 

„, 2GJcosi. 2GJcosi (1 + e cos f ) . 
W K gm = --r—u = — \ {l _ e ^ u. 

In eq.(27) u ~ / is assumed due to the the fact that the osculating element ui may be 
retained almost constant on the temporal scale of variation of the true anomaly of the 
test body. 



4. Secular gravitomagnetic effects on the Keplerian orbital elements: spher- 
ical central source 

The secular effects can be worked out by adopting the same strategy followed in [16] 
for a similar kind of perturbing functions. When eq.(27) is averaged over one orbital 
period P of the test body, a, e, i, f2 and w are to be considered constant: 

„ 1 [ P G2J , n GJcosi f 27T df n G Jcosi 

The relation du = duo + df = df has been used in eq.(28) which can now be used in 
determining the secular changes of the Keplerian orbital elements of the test body. It 
yields: 



(29) 


d< n > 2n 


dM 


(30) 


d< 1Z >2-k 


duj 


(31) 




on 


(32) 


d< TZ > 2 K 


di 


(33) 




de 



^ G J sin i 

c 2 a(l — e 2 ) ' 

G Je cos i 

— — 4ti 

c 2 a(l — e 2 ) 2 
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A particular care is needed for the treatment of n when the derivative of < 1Z >i- K with 
respect to a is taken; indeed, it must be posed as: 



(34) 

In eq.(34) 

(35) 
and 
(36) 



d< K >2n d< TZ > 2 K d< TZ >2n dn 



da 



da 



In 



+ 



dn \ a da' 



d<lZ>2Tr „ G Jcosi 
= 2n 



da 



c 2 a 2 (l 



2V 



d<lZ>2TT dn G Jcosi 
dn \ a da c 2 a 2 (l — e 2 ) 



From eq.(17) and eq.(29) it appears that there are no secular changes in the semimajor 
axis, and so the orbital period of the test body, related to the mean motion by P — 2ir/n, 
can be considered constant. This implies that in eq.(34) only eq.(35) must be retained. 
Using eqs.(29)-(34) in eqs.(17)-(22) one obtains for the secular rates: 



(37) 
(38) 
(39) 
(40) 
(41) 
(42) 



da 
d* 
de 
di 
di_ 
dt 

dn 
~dt 

dui 
~dt 
dM 
dt 



LT 



LT 



LT 



LT 



LT 



LT 



o, 
o, 
o, 

G 



2J 



c 2 a 3 (l 
G 



. e 2)3/2 - 

6J 



c 2 a 3 (l - e 2 f/ 2 



cost, 



0. 



They are the well known Lense-Thirring equations [1, 2]. In deriving them it has been 
assumed that the spatial average over / yields the same results for the temporal average 
[17]. 



5. Secular gravitomagnetic effects on the Keplerian orbital elements: non- 
spherical central source 

In this Section we shall deal with a non spherical central rotating source with axial 
symmetry around the rotation axis. 

In classical electrodynamics the potential vector for a generic steady current distri- 
bution can be written as: 



(43) 



A(r) = i/ 1 ^>dr'. 
c J | r-r | 
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The quantity p is the charge density which, in general, depends on r, but, in this case, 
not on time. If we assume that the current distribution rotates uniformly around an axis, 
chosen as z axis, then, since for any current element v = az x r with a angular velocity 
of the current distribution eq.(43) becomes: 

(44) A(r) = -J j^-^dr . 

In order to deal with an axisymmetric current distribution let us introduce the cylindrical 
coordinates: 



(45) 



£ > 0, < < 2n. In this case it can be shown that the lines of the potential vector are 
circles around the z axis and its modulus is a cylindrical symmetric function of £ and z. 
So, eq.(44) reduces to: 

(46) A(r) = -H(£, z)e 4> . 

The function if(£, z) remains unchanged under rotations around the z axis. Passing 
from cylindrical to rectangular Cartesian coordinates, the components of the potential 
vector become: 

(47) A l = --H(£, z)y, 

c 

(48) A 2 = -H{£, z)x, 

c 

(49) = 0. 

Recalling that, in the linear approximation, the general relativistic gravitomagnetic po- 
tential can be obtained from the vector potential of electromagnetism times —AG [2], 
eqs.(47)-(49) lead to: 

(50) h 01 = —H{£, z)y, 

4Ga 

(51) h Q2 = z)x, 

c 

(52) h 03 = 0. 

Of course, eqs.(50)-(52) can be rigorously obtained solving the linearized Einstein field 
equation, written in the Lorenz gauge, for a localized (g^ u — > i]^ v at spatial infinity), 
axisymmetric, stationary mass distribution in the weak field and slow motion approxi- 
mation: 

(53) A^ ofc =0, k= 1,2,3 

. i . , 167rG m , „ „ 

(54) Ah ok - -^T ok , k = 1,2,3. 
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The quantities T fc are the {Ofc} components of the stress-energy tensor for the matter; 
in deriving eq.(54) the internal stresses have been neglected. Eq.(54) is valid inside the 
matter, while eq.(53) holds in the free space outside the central body and tells us that 
H(£, z)y and H(£, z)x are harmonic functions. This feature was used by Teyssandier 
in obtaining a multipolar expansion of H(£, z) [11]. Introducing in the frame K the 
usual spherical coordinates {r, 9, <f>}, if r e is the radius of the smallest sphere centered 
on the origin of the coordinates containing the whole body (in practice, it should be the 
equatorial radius, _R® in the case of the Earth), in the region r > n > r e H(r, 9) is 
given by: 

j oo 

(55) H(r, e) = —[l-J2K l (-) l Pi + i(cos9)l 

zr i=i r 

with: 

2 Mr 2 

(56) ^=27T3^ (W <+ 2) ' 



(57) U = / P( r '> ^>' Z+2 fl(cos0')dr'. 

In eqs.(55)-(57) / is the moment of inertia of the body about the z axis, M is the total 
mass of the central body, p is its density, P J+1 (cos 9) is the first derivative of the Legendrc 
polynomial of degree I + 1 and Ji is the Newtonian multipole moment of degree I, given 
by: 

(58) J l = -A-J p (r, 9')r' l p l (cos9')dr'. 

It is interesting to note that, if in eq.(55) only the spherical symmetric term: 

(59) H^(r) = ± 

is retained, as it would be the case for a perfectly spherical central body, eqs.(50)-(52) 
reduce to: 

(60) <» - 

<«') C - ~ 

(62) h'S = 0. 



For a spherical rotating body J = la, and so eqs.(60)-(62) can be cast into the familiar 
form: 

m a<°> . -H£4i. 
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The correction of order I to the gravitomagnetic potential due to the nonsphericity of 
the central rotating body is given by: 

(64) f$ = y K^yP^icosO), 

(65) f $ = ?^ x K l Cf) l Pl +1 (co80), 

(66) h$ = °- 

Modeling the central body as a spheroid stratified into ellipsoidal shells, in [11] is shown 
that all the relativistic coefficients Ki are null, except K2', if we consider our planet, from 
the analysis of Teyssandier of various models of Earth's interior, it results to be positive 
and of 0(1CT 3 ). 

The starting point in deriving the relativistic multipole corrections of degree I to 
eqs.(37)-(42) is the perturbative lagrangian term: 

(67) C% = ™(A(0.v), 

where m is the mass of the point particle and v is its velocity. In the case / = 2, in 
eqs.(64)-(66) the first derivative of the Legendre polynomial of degree 3 appear; 

(68) P 3 (q) = l - (5q 3 -3q), q^cosO, 

(69) P s {q) = l - (15<z 2 -3), q^cosO. 

Inserting eq.(69) in eqs.(64)-(66) for I = 2 leads to: 

(70) Cfi = -m°g rlK2{15co8 i _ 3). 

Such a perturbative term must be expressed in term of the Keplerian orbital elements in 
order to be employed in the first order Lagrange planetary equations. In this case TZ is 
TZ^ = ■ v). Using eq.(15) and recalling that z — rcos9 it is possible to obtain: 

(71) = — ucosi 

(72) cos 8 = sin u sin i. 

Neglecting terms of order 0(e n ), n > 2, we can write: 

1 ~ 1 + 3ecos / 

{ ' r 3 ~ a 3 (l - e 2 ) 3 ' 

having assumed: 

1 1 + e cos / 



(74) 



r a(l — e 2 ) 
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Eqs.(71)-(73) in cq.(70) give: 
(75) 



_ r2 i Gla r 2 K 2 cosi . 2 . 2 . 

TZ 1 - 2 ' = — —-In — (1 + 3ecos/)(15sm 2 w sm 2 i - 3)u. 



c 2 a 3 (l-e 2 ) 3 

If we want to investigate the secular trends of the orbital elements, we must average 
eq.(75) over an orbital period of the test body, as done in the previous Section. It is 
straightforward to obtain: 

(76) < (l + 3ecos/)(15sin 2 wsin 2 i-3) u > 27r = |(9 - 15cos 2 i). 

In deriving eq. (76) we have adopted the reasonable assumption that the pericenter of the 
test body remains almost unchanged during an orbital revolution, i. c. du = duo+df ~ df. 
So we have: 



(77) 



< K (2) >2n= n 



Gla r 2 e K 2 
2c 2 a 3 (l-e 2 ) 3 



(9cosi — 15cos 3 i). 



Eq. (77) can be considered as the 2nd order correction to the gravitomagnetic perturbing 
function given in eq.(28). Eq.(77) yields: 



(78) 
(79) 
(80) 
(81) 
(82) 
(83) 



d< 


>2tt 


dM 




d< 


>2tt 


du 






>2tt 


an 




d< Tl^ 


>2it 


di 




d< -rw 


>2tt 


de 



d< > 27T 



= 0, 
= 0, 
= 0, 
= n- 



Gla r 2 e K 2 
2c 2 a 3 (l-e 2 ) 3 
ieGIa r 2 K 2 



sini(45cos 2 i — 9), 



c 2 a 3 (l-e 2 ) 4 
SGI a r 2 K 2 



da 



I, 



2c 2 a 4 (l-e 



2y;> 



(9cosz — 15cos 3 i), 



(9cosz — 15cos 3 i). 



By inserting eqs.(78)-(83) in eqs.(17)-(22) it can be obtained for the secular rates: 
(84) 



(85) 
(86) 
(87) 



da 


( 2 ) 


dt 


LT 


de 


( 2 ) 


dt 


LT 


di 


( 2 ) 


dt 


LT 


dn 


( 2 ) 


~dl 


LT 



2c 2 a5( 1 „ e 2)7/2 



(45cos 2 i - 9), 
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cLv 
d* 



(2) 



. dn w 3GIa 





(88) 



= — COS I — — 



c 2 a 5 (l - e 2 ) 7 /2 



(9cosi — 15cos 3 i), 



LT 



AM 



(2) 



(89) 



= 0. 



LT 



Eqs.(84)-(89) yield the corrections to the precessional rate Lense-Thirring equations when 
the central rotating body is not perfectly spherical but only axially symmetric. 

If we use in eqs.(87)-(88) the values r e = R e ~ 6378 km, K 2 = 0.874 x 10~ 3 [12] we 
can obtain an estimate of the sensitivity of LAGEOS and LAGEOS II to the relativistic 
Earth's quadrupole correction to the Lense-Thirring precessional rates. The results are: 

(90) ^lageos = ~6.7 x 10 _1 mas/century, 

(91) ^lageosii = 2 - 8 mas/century, 

(92) ^lageosii = 5 - 4 mas/century. 

The present accuracy in the measurement of the LAGEOSs' rates of the node and the 
perigee is of the order of 1 mas/y. This implies that such tiny corrections do not affect 
the current efforts in detecting the Lense-Thirring drag. 

6. Conclusions 

In the weak field and slow motion approximation of general relativity we have de- 
veloped an alternative approach to the calculation of the Lense-Thirring effect on the 
Keplcrian orbital elements of a test particle freely falling in the field of different kinds of 
axially symmetric central sources. 

Such strategy stresses the formal analogy with the electromagnetism recalling the 
lagrangian of a charged particle acted upon by the Lorentz force and create a link to 
space geodesy by exploiting the widely used Lagrangian planetary equations. 

For a perfectly spherical source the well known Lense-Thirring processional rates for 
the node and the perigee are obtained. 

If departures from sphericity of the gravitational source are accounted for, it has been 
found that only the node and the perigee are affected by them through additional secular 
rates. 

At present they cannot influence the current measurement of the Lense-Thirring drag 
by means of LAGEOS and LAGEOS II because they fall below the experimental sensi- 
tivity. 



I am profoundly indebted to I. Ciufolini who focused my attention to the fascinating 
topic of gravitomagnetism. I would like also to thank M. Gasperini for discussions, P. 
Colangelo who has read carefully the manuscript and L. Guerriero who supported me at 
Bari. 
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